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Abstract. Shape-dependent universal crossing probabilities are studied, via Monte Carlo simulations, for bond and 
site directed percolation on the square lattice in the diagonal direction, at the percolation threshold. In a dynamical 
interpretation, the crossing probability is the probability that, on a system with size L, an epidemic spreading without 
immunization remains active at time t. Since the system is strongly anisotropic, the shape dependence in space-time 
enters through the effective aspect ratio r^g — ct/L^, where c is a non-universal constant and z the anisotropy exponent. 
A particular attention is paid to the influence of the initial state on the universal behaviour of the crossing probability. 
Using anisotropic finite-size scaling and generalizing a simple argument given by Aizenman for isotropic percolation, 
we also obtain the behaviour of the probability to find n incipient spanning clusters on a finite system at time t. The 
numerical results are in good agreement with the conjecture. 

PACS. 64.60.Ak Renormalization-group, fractal, and percolation studies of phase transitions - 05.50.-l-q Lattice 
theory and statistics (Ising, Potts, etc.) - 02.50.-r Probability theory, stochastic processes, and statistics 



1 Introduction 

. Most of the properties of the critical state at a second order 
' phase transition are related to the existence of self-similar cor- 
related clusters. In the percolation problem, the critical point is 
characterized by the spontaneous occurence of infinite clusters 
of connected sites, crossing the system in the thermodynamic 
. limit. 

The study of crossing probabilities for standard percolation 
, on finite systems at the percolation threshold has been an ac- 
' tive field of research during the last decade [1-44] (for a recent 
' review see |45 1). 

In two dimensions, the crossing probability tt can be de- 
fined as the probability to find at least one cluster joining two 
opposite edges of a rectangular-shaped system with length i|| 
and width L±. At the percolation threshold, in the finite-size 
scaling (FSS) limit, i.e., when L|| oo, L± oo, while the 
. aspect ratio r — L^^/L± remains fixed, the crossing probabil- 
ity is a scale -invariant universal function ^{r) of the aspect 
ratio jm. 

Following the extensive numerical study of Langlands et 
al. ID, an analytical expression for 7r(r) at criticality was de- 
rived by Cardy |2 1. Using the relation between percolation and 
the cluster formulation of the q-state Potts model in the limit 
q ^ 1 [46 '47 1 as well as the techniques of boundary confor- 
mal field theory, he was able to calculate the probability for 
two non-overlapping segments on the edge of the half-plane to 
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be connected. The crossing probability in the rectangular ge- 
ometry could then be deduced from the appropriate conformal 
mapping. One may notice that the scale invariance of 7r(r) is a 
non-trivial property, resulting from the occurence of a vanish- 
ing scaling dimension x{q) associated with a boundary condi- 
tion changing operator of the Potts model in the limit q ^ 1 
for percolation. 

Other exact crossing formulas were later obtained for dif- 
ferent geometries |5 6|. Since then, mathematically rigorous 
proofs have been provided for some of these results [25 26l. 

Crossing probabilities on same-spin Ising clusters [27 28 29| 
and Potts clusters L30J in two dimensions have been also con- 
sidered. 

Recently, some attention has been paid to the properties of 
incipient spanning clusters [31-44]. The probability to find n 
incipient spanning clusters, connecting two disjoint segments 
at the boundary of a finite critical system, has been studied 
through Monte Carlo simulations [31-39]. Rigourous bounds 
on the spanning probability |40 41 1 have been obtained as well 
as analytical expressions through conformal and Coulomb-gas 
methods II42I43I44I . 

In two recent works f48"491, we started a study of the criti- 
cal crossing probability (in space-time) for off-equilibrium sys- 
tems in 1h-1 dimensions. In these strongly anisotropic systems, 
when the temperature-like scaling field A is non-vanishing, the 
correlation functions generally display an exponential decay, 
in the space and time directions, characterized by a correlation 
length ^ and a relaxation time r diverging respectively as I.50J 

C = eZ\-% r = TA-^\ (1.1) 
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when A ^ 0. The prefactors ^ and f are non-universal ampli- 
tudes and z is the dynamical (anisotropy) exponent. 

The directed percolation problem |51 52 53| was studied 
numerically through Monte Carlo simulations in |48 1 and ana- 
lytical results were obtained for an exactly solvable diffusion- 
limited coagulation process in 1491 . In these sytems, the cross- 
ing probability in the time direction vTf is defined as the proba- 
bility that the system of size L remains active at time t. Aniso- 
tropic scaling 1501541 then implies that the appropriate aspect 
ratio is r = t/L'', i.e., the rescaled time. We found that in the 
FSS limit, as for isotropic critical systems, the critical crossing 
probability is a scale-invariant universal function of an ejfec- 
tive aspect ratio which is the product of r by a non-universal 
constant c l55l . 

In the present work, we continue the examination of the 
universal properties of the crossing probability ttj in directed 
percolation which, in a dynamical interpretation, essentially 
corresponds to the survival probability of an epidemic spread- 
ing without immunization |56|. Sectionals devoted to a de- 
tailed study of the influence of the initial state on the universal 
behaviour of ttj . In Section|3]we present some results concern- 
ing the scaling behaviour of the probability to find n incipient 
spanning clusters, i.e., n disjoint centres of infection surviving 
in the FSS hmit. 




x=l 2 L-l L 



Fig. 1. Directed percolation on the square lattice in the diagonal direc- 
tion. There are L sites in the spatial direction which are either occu- 
pied (•) or empty (o). Three types of initial states are considered with 
a) L occupied sites, b) / < L nearest-neighbour occupied sites and 
c) a single occupied site. In the last two cases, for FBC, the occupied 
sites are located in the middle of the system. 

7rt(r, /), in the FSS limit where L, t and I — * oo, while keeping 
r and / fixed. 



2 Critical crossing probability 

We study the critical crossing probability Pt, in the time di- 
rection, for bond and site directed percolation in 1 + 1 dimen- 
sions through Monte Carlo simulations on the square lattice. 
The time axis is oriented in the diagonal direction. The space 
and time coordinates take alternatively integer and half-integer 
values on successive spatial rows. On a site located at [x, t), 
with \ < X < L, two directed bonds are leaving, which ter- 
minate on the nearest-neighbour sites at [x ± 1/2, t + 1/2) 
as shown in Figure [J We use either free boundary conditions 
(FBC) or periodic boundary conditions (PBC) in the spatial di- 
rection. In the later case, x is defined modulo L. 

In the epidemic language, occupied (empty) sites corre- 
spond to infected (healthy) individuals. In the bond problem, 
the bonds are independently open with probability p. A site is 
occupied at t + 1/2 when it is connected via an open bond to 
a site which was occupied at time t. In the site problem all the 
bonds are open and a site is occupied with probability p at time 
f + 1/2 when at least one of its first neighbours at time t is 
occupied. 

A directed percolation cluster is a collection of connected 
occupied sites starting from some source at i = 0. We con- 
sider successively the three types of initial states shown in Fig- 
ure [2 A sample contributes to the crossing probability when 
at least one cluster survives at time t. The simulations are per- 
formed at the percolation threshold, pj?""'^ = 0.644700185(5) 
for the bond problem and = 0.70548522(4) for the site 
problem. These values, as well as the dynamical exponent z = 
1.580745(10), are taken from |57|. The crossing probability 
Pt(L, t, I) is determined as a function of the aspect ratio, r = 
t/L^, and the fraction / = l/L of occupied sites in the ini- 
tial state, using 10^ samples with sizes up to L = 2^". When 
/ is non-vanishing, Pt converges to a scale-invariant function. 



2.1 Initial state with a fully occupied lattice 

When the L sites are randomly occupied with probability pi > 
0, the crossing probability converges to a scale-invariant func- 
tion in the FSS limit. This is illustrated in Figure|2lin the case 
of site percolation with FBC. When pi — 1 the convergence 
is from above and the results obtained with L = 8 are already 
quite close to the FSS limit. Thus pi is an irrelevant variable 
when / > 0. In the following we always take pi — 1 in or- 
der to obtain a better convergence towards the scale-invariant 
behaviour 

We now consider site and bond directed percolation with 
either FBC or PBC in order to study the universality of the 
crossing probability. The lattice is fully occupied in the initial 
state. The behaviour of the crossing probability as a fonction 
of the aspect ratio is shown in Figure |3] The raw data con- 
verge quickly to scale-invariant functions as the system size 
is increased. They are shown in the inset for the largest size, 
L = 256. Different functions are obtained for the two types of 
boundary conditions and the curves for the bond and site prob- 
lems are shifted with respect to one another, along the horizon- 
tal axis, by an amount which is the same for FBC and PBC. 
The main figure shows the data collapse on universal curves, 
one for each type of boundary condition, after the results for 
site percolation have been shifted. 

In order to explain this behaviour let us consider the scaling 
behaviour of the crossing probability for an off-critical system 
with a deviation A = \p ~ Pc\ from the percolation threshold. 
Under an anisotropic change of the length scale by a factor b 
the length transforms as i' = L/b, the time ast' ^ t/b^ and A 
as A' = b^/^A, where 1/ is the coiTelation length exponent de- 
fined in Assuming that, as indicated by the Monte Carlo 
data, the crossing probability is scale invariant, then it satisfies 
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Fig. 2. Crossing probability as a function of In r for site percolation 
with FBC. With increasing size, the values of nt obtained when the L 
sites are randomly occupied in the initial state with probability pi = 
1 /4, converge to the values obtained for = 1 on a lattice with size 
L = 256. The statistical errors are smaller than the symbols. 



the scaling relation 

Pt{L,t,A)=Pt 

With b = A^'^ one obtains 
L 



Pt{L,t,A) = Pt 



A-^' A- 



1 = 



L t 



(2.1) 



(2.2) 



where g{u, v) is a universal fonction of its dimensionless argu- 
ments. The FSS limit corresponds to Z\ = and b — L\n \2.\\ . 
which gives 

Pt(L,i,0) -Ft (^1,-^,0^ -7rt(cr). (2.3) 

For each type of boundary conditions, the crossing probability 
is a scale-invariant universal function of the effective aspect 
ratio Tcff — cr. The non-universal amplitude c depends on the 
choice of the length and time units. It can be expressed as a 
function of ^ and f |54|. Comparing (|^} and ( 12. 3> , r appears 
through the dimensionless ratio v ju^ so that 



(2.4) 



The non-universal constant c is different for the two percolation 
problems. According to ( 12. 3> . for a given type of boundary con- 
dition, the crossing probabilities TTt are identical for the site and 
bond problems when the values of the aspect ratio, and rb, 
satisfy the relation c^r^ ~ Chfh- Thus the data collapse in Fig- 
ure|3]is obtained through a shift of the site percolation data by 



In Tb — In Ts 



In 



Cb 



.2476(8) 



(2.5) 



The numerical value of S^^ was estimated in |48| through a 
least-square fit. 



Fig. 3. Crossing probability as a function of In r with a fully occu- 
pied lattice in the initial state. The Monte Carlo results for the largest 
size (L = 256) are shown in the inset for the bond and site problems 
with FBC and PBC. The main figure shows the data collapse on two 
different universal curves for FBC and PBC in the scaling limit, af- 
ter shifting the data for site percolation, as explained in the text. The 
statistical errors are smaller than the symbols. 



The crossing probability can be calculated using a transfer 
matrix in the time direction T, working with a restricted basis 
of states corresponding to configurations with at least one oc- 
cupied site |58 1. The matrix element Taji gives the probability 
that the configuration in state |a) at time t leads to the config- 
uration in state |/9) at time t + 1. Given an initial state the 
crossing probability then reads 



(2.6) 



When r — t/L^ ^ 1, the leading contribution comes from the 
largest eigenvalue ylmax < 1 of T and 



exp 



ct 



exp(— acr) 



(2.7) 



since the relaxation time | ln(yl„iax)|^^ scales as L^. In the 
argument of the exponential, a is some universal constant. Thus 
InTTf is a universal linear function of cr when r 3> 1. This 
universal linear variation is shown in Figurel^where r for the 
site problem has been multiplied by Cg/cb = 1.281(1) in order 
to eliminate the influence of the non-universal constant c on 
the data collapse. Since the argument leading to ( 12. 7> is quite 
general, the same behaviour is expected for the other types of 
initial states considered below. 



2.2 Initial state with a sequence of occupied sites 

Let us consider now the behaviour of the crossing probabil- 
ity when I consecutive sites are occupied in the initial state, 
as shown in Figure^. The results for PBC and FBC are pre- 
sented in Figures |5l and |6j respectively, for three values of the 
occupation ratio / = l/L, from top to bottom / ~ 1/2, 1/4, 
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Fig. 4. Logarithm of the crossing probability as a function of r with 
a fully occupied lattice in the initial state. For the site problem, r has 
been multiplied by Cs/c]^. In the scaling limit, a linear universal be- 
haviour is obtained, with different slopes for FBC and PBC. 

1/16. The raw data for the largest size, L ~ 256, are shown in 
the insets. As above, different curves are obtained for the site 
and bond problems. Shifting the site results by J^j,, the data 
collapse on three universal curves, one for each value of /, as 
shown on the main figures. 

The problem involves the new length scale, I, and the cross- 
ing probability scales as 



Pt{L,t,l,A)^Pt 

With b = A^'^ one obtains 
L 



L LI b^-A 



(2.8) 



PtiL,t,l,A)=Pt 



(2.9) 

where u, w) is a universal function of its arguments. In the 
FSS limit, at the critical point, equation (12. 8> gives 

Pt{L,t,l,0)^Pt (^1,-1,1,0^ =7Tt{crJ). (2.10) 

The dimensionless occupation ratio / does not introduce 
any new non-universal constant in the scaling limit where L 
and I oo. Thus, as shown in Figures |5] and |6l the crossing 
probability is a scale-invariant universal function of the effec- 
tive aspect ratio Vcff, which depends on the value of / as well as 
on the boundary conditions. One may notice a slightly slower 
convergence to the limiting behaviour when / decreases. 



2.3 Initial state with a single occupied site 

When a single site (located in the middle of the system for 
FBC) is occupied in the initial state, the crossing probability 
TTj^-* is no longer scale invariant. It acquires a dimension which 
is the scaling dimension of the order parameter. 



V 



.252072(8), 



(2.11) 



Fig. 5. Crossing probability as a function of In r with a fraction / = 
l/L of occupied lattice in the initial state and PBC. The Monte Carlo 
results for the largest size (L = 256) are shown in the inset for the 
bond and site problems with / = 1/2, 1/4 and 1/16 from top to 
bottom. The main figure shows the data collapse in the scaling limit 
on three different universal curves, one for each value of /, when the 
data for site percolation are appropriately shifted. The statistical errors 
are smaller than the symbols. 




Fig. 6. As in Figure|5]for FBC. 

corresponding to the values 

13 = .276486(8) , v = 1.096854(4) , (2.12) 

taken from f57V This behaviour follows from the scaling of 

p/^\ the probability that a cluster grown from a single central 
occupied site survives after t time steps. 



(1) 



L t 



b^/^A 



(2.13) 



On the infinite system with b = t^/^, one obtains 



Pi^^\oo,t,A) = i-^-/^pf ) (oo,l,ti/^'^Z\) , (2.14) 



in agreement with equation (94) in reference 1531 . where the 
exponent (3' /t^w = (3/v\\ = Xm/ z for directed percolation with 
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Fig. 7. Crossing probability for site (dotted line) and bond (full line) 
directed percolation with a single occupied site in the initial state as 
a function of Inr. After rescaling tTj^' by a factor L^™, the size ef- 
fects are eliminated. Although the data for the site problem have been 
shifted by S^^y^, a different behaviour is still obtained for site and bond 
percolation. 
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Fig. 8. Once ln(L^'" vr^^' ) for the site problem has been also shifted by 
5sb in order to eliminate the influence of the non-universal amplitude 
in <2.19> . the data collapse on two universal curves for FBC and PBC, 
respectively. For small values of In r a linear variation is obtained with 
a slope —Xm/z, as indicated by the dashed line. 



our notations. Taking b = L at criticality (A — 0) leads in the 
FSS limit to 



(2.15) 



the crossing cluster issued from the middle of the sample does 
not explore the boundary regions. For the same reason Trj^-* is 
independent of the size L in this regime. Different curves are 
obtained for the site and bond problems for a given type of 
boundary conditions, meaning that a new non-universal con- 
stant is involved. In a log-log plot, it is clear that a single curve 
can be obtained through a vertical shift (5^^^ of the data for the 
site problem. The resulting data collapse is shown in Figure|8] 
Thus the new non-universal constant A is a multiplicative one, 
either As or A\,, and the shift is given by: 



^^b - In 



.178(4) 



(2.16) 



It was calculated by first shifting the data of the site problem 
by (5g^^. Due to this horizontal shift of the site data, a direct cal- 
culation of the vertical shift is not possible. Thus, considering 
y = ln(L^™7r(^^) as a function Inr, the data of the bond prob- 
lem were fitted to a quadratic polynomial in the linear region, 
for small values of In r. Then S^-^^ was obtained as the average 
of Vh — Us, evaluated at small enough values of In r for the site 
problem. 



Alternatively, the scaling behaviour of tt^ 



I can be 



recovered by assuming that, asymptotically when / <C 1, the 
scale-invariant crossing probability in (I2.10> behaves as 



(2.17) 



lim TTt (cr, /) = ip{cr) 

so that L^^TTf ~ F"^ip{cr). 

Furthemore the non-universal amplitude, Ag or Ab, which 

appears only in ttj^'*, i.e., when the number I of occupied sites 
in the initial state is 0(1), can be explained if one associates to 
these / sites an effective length l^g given by 



'cff = I + 61 , 61 — 6ls, 6lh 



(2.18) 



where the increment 61 is 0(1) and non-universal. Thus, replac- 
ing / by /cff = {I + 61) /L, there is no change in the FSS limit 
when / > 0, but when / = 1 we obtain 



i^-TT^^^ = {l + 6lf'^<fi{cr) 
with a non-universal amplitude given by 

A={l + 6lf"^ . 



(2.19) 



(2.20) 



As shown in Figure Q a good data collapse is obtained for numerical results in FigurelH 



According to (I2.16> we have 6ls > 6lh. 

Since, as noticed above, tt^^'' becomes independent of L 
for small values of r, the scaling function (p{w) with w = 
cr has to behave as w"'^"^!^ when w <^ \. This leads to a 
slope —Xm/z = —fi/zv — —.159464(5) for the linear part of 

ln(L^™7rj^^) as a function of Inr, in good agreement with the 



In 



2^£c„^(i) ^jjjj jjjg different system sizes, L = 2'^ to 2 
this figure, the data for the site percolation problem have been 
shifted by 6sh given in \2.5\ in order to take into account the 
non-universal constants Cs and Cb. 

For small values of r, » t, the boundary conditions are 
irrelevant since, with a single occupied site in the initial state. 



3 Critical spanning clusters 

In this Section we present a preliminary study of the scaling be- 
haviour of the probability 7rt„ (r) to find exactly n independent 
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Fig. 9. Probability ntn to find n critical spanning clusters in the time 
direction for the site problem with FBC when L = 256. The number 
of clusters is n = 1 to 4 from right to left. For n = 1, the spanning 
probabilities obtained when L = 64 and 128 are also shown. 



clusters surviving after t = rL^ time steps on a system with 
size L in the FSS limit. In the epidemic interpretation iTtn is 
the probability to find n disjoint centres of infection surviving 
in the FSS limit. This problem presents some similarity with 
the statistics of family names first studied by Gallon and Wat- 
son t59J (see also |60 61|). 

The spanning probability was estimated by generating 10^ 
samples with L = 64, 128 and 256 for the site problem with 
FBC or PBC and L occupied sites in the initial state. As the 
size increases, it converges to a scale-invariant function 7rt„ of 
the effective aspect ratio r^s- A better convergence is obtained 
with FBC, thus we shall only analyse the results obtained with 
this type of boundary conditions. An example of the raw data 
for TTtn at the largest sizes studied is shown in Figure|9l 

The scaling behaviour can be obtained by extending to the 
case of a strongly anisotropic system a simple argument pro- 
posed earlier for isotropic percolation [40 4 2 |. 

Let Pt (n, r) denote the probability to have exactly n span- 
ning clusters at time t, on a system in c? + 1 dimensions, with 
volume L'^. The scale invariance of Pt implies that the size de- 
pendence is wholly contained in the aspect ratio r. Dividing 
the system in subsystems with the same volume (L/b)^, as 
shown schematically in Figure^] for d = 1 and b = n, the 
dominant event will correspond to rt' = n/b'^ spanning clus- 
ters in space-time for each subsytem. Since the aspect ratio of 
the large system is r = i/L^, it becomes r' = t/{L/by = rV- 
in the smaller subsystems. Thus, up to prefactors, one expects 
the following behaviour: 



(3.1) 



According to the transfer-matrix argument of Section 12711 a 
leading exponential dependence on t (through r) is expected. 
Thus, with b — ■n}/'^, we obtain 





Ub 



Fig. 10. n spanning clusters on a system with size L at time t in 1 + 1 
dimensions are considered to result typically from the juxtaposition of 
n/b spanning clusters on b subsystems with size L/b. 



Applying the first relation to the exponential form allows us to 
identify the exponent uj = l + z/dso that the scale-invariant 
crossing probability behaves as 



■Ktn{r) ~ exp(-arn^+^/'') 



(3.3) 



Pt{n,r)^ Pt(l,W'*r) ^ exp( 







(3.2) 



where a is a non-universal prefactor Actually, the crossing pro- 
bability TTtn, like TTj, is a universal function of the effective as- 
pect ratio cr. This scaling behaviour generalizes the form con- 
jectured by Aizenman for isotropic percolation m D = d + 1 
dimensions |40J, which was recently tested numerically with 
success |39|. The isotropic limit corresponds to z = 1 and 
d = D — I 'm ( 13. 3> . The anisotropic result should remain valid 
below the upper critical dimension c?c = 4 for directed percola- 
tion 1 52 1. At and above d^, as in the isotropic case, one expects 
a proliferation of spanning clusters, linked with the breakdown 
of hyperscaling |4()"62l. 

The Monte Carlo results in Figures ^3 and El for the site 
and bond percolation problems with FBC at the largest size, 
L — 256, are in good agreement with this scaling behaviour. 
The insets show the linear variation of In iitn when r is suf- 
ficiently large. Different slopes are obtained for different val- 
ues of n. As shown in the main figures, after multiplication by 
good data collapse is obtained. 

The universality of iTtn [fcs) is illustrated in Figure[21 The 
dotted line on the main figure corresponds to the site problem 
with n = 2 after rescaling r by the non-universal factor Cg/cb. 
The remaining shift between the site and bond data is likely to 
be due to finite-size corrections. 



4 Conclusion 

For the directed percolation problem in the FSS Umit, the cross- 
ing probability in the time direction, ttj, is a scale-invariant uni- 
versal function of the effective aspect ratio r^s — ctjU' , ap- 
propriate for a strongly anisotropic system. It also depends on 
the fraction f — 1/ Lof occupied sites in the initial state, where 
I is the length of the sequence of occupied sites. This function 
vanishes with / as /^™, where Xm — (3/i^ is the scaling di- 
mension of the order parameter. When the L sites are occupied 
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cient simulation method, like the "go with the winner" strategy 
introduced by Grassberger r36"38"391. It would be particularly 
interesting to study the expected crossover of TTm, due to the 
proliferation of clusters, at and above the upper critical dimen- 
sion 1 39 1. It should be easier to do this for directed percolation, 
where dc = 4, than for isotropic percolation, where Dc = 6. 

The present work can be also extended by considering the 
crossing probability on same-spin clusters in strongly anisotro- 
pic spin systems at equilibrium, like the three-dimensional uni- 
axial ANNNI model at its Lifshitz point, for which the critical 
parameters are now accurately known L63J . 
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Fig. 11. Scaling behaviour of the crossing probability ntn(r) for the 
site percolation problem with FBC when L = 256. The inset shows 
the linear variation of In -Ktn versus r with a slope depending on n = 2 
to 6 from top to bottom. A good data collapse is obtained in the main 
figure, in agreement with <3.3> for d — 1, when In ixtn is rescaled by 




Fig. 12. As in Figure fTTI for the bond percolation problem. The dotted 
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